Microbial Dynamics and Pulmonary Immune Responses in COVID-19 Secondary Bacterial Pneumonia (preprint)

Secondary bacterial pneumonia (2°BP) is associated with significant morbidity following respiratory viral infection, yet mechanistically remains incompletely understood. In a prospective cohort of 112 critically ill adults intubated for COVID-19, we comparatively assessed longitudinal airway microbiome dynamics and studied the pulmonary transcriptome of patients who developed 2°BP versus controls who did not. We found that 2°BP was significantly associated with both mortality and corticosteroid treatment. The pulmonary microbiome in 2°BP was characterized by increased bacterial RNA load, dominance of culture-confirmed pathogens, and lower alpha diversity. Bacterial pathogens were detectable days prior to 2°BP clinical diagnosis, and in most cases were also present in nasal swabs. Pathogen antimicrobial resistance genes were also detectable in both the lower airway and nasal samples, and in some cases were identified prior to 2°BP clinical diagnosis. Assessment of the pulmonary transcriptome revealed suppressed TNFa signaling via NF-kB in patients who developed 2°BP, and a sub-analysis suggested that this finding was mediated by corticosteroid treatment. Within the 2°BP group, we observed a striking inverse correlation between innate and adaptive immune gene expression and bacterial RNA load. Together, our findings provide fresh insights into the microbial dynamics and host immune features of COVID-19-associated 2°BP.

Comprehensive characterization of the transcriptional response to COVID-19 in multiple organs reveals shared signatures across tissues (preprint)

Infection by Coronavirus SARS-CoV2 is a severe and often deadly disease that has implications for the respiratory system and multiple organs across the human body. While the effects in the lung have been extensively studied, less is known about COVID-19 cellular impact across other organs. Here we contribute a single-nuclei RNA sequencing atlas comprising six human organs across 20 autopsies where we analyzed the transcriptional changes due to COVID-19 in multiple cell types. Computational cross-organ analysis for endothelial cells and macrophages identified systemic transcriptional changes in these cell types in COVID-19 samples. In addition, analysis of signaling pathways from multiple datasets showed several systemic dysregulations of signaling interaction in different cell types. Altogether, the COVID tissue atlas enables the investigation of both cell type-specific and cross-organ transcriptional responses to COVID-19, providing insights into the molecular networks affected by the disease and highlighting novel potential targets for therapies and drug development.

Robust Variation in Infant Gut Microbiome Assembly Across a Spectrum of Lifestyles (preprint)

Infant microbiome assembly is intensely studied in infants from industrialized nations, but little is known about this process in populations living non-industrialized lifestyles. In this study we deeply sequenced infant stool samples from the Hadza hunter-gatherers of Tanzania and analyzed them in a global meta-analysis. Infant microbiomes develop along lifestyle-associated trajectories, with over twenty percent of genomes detected in the Hadza infant gut representing phylogenetically diverse novel species. Industrialized infants, even those who are breastfed, have microbiomes characterized by a paucity of Bifidobacterium infantis and gene cassettes involved in human milk utilization. Strains within lifestyle-associated taxonomic groups are shared between mother-infant dyads, consistent with early-life inheritance of lifestyle-shaped microbiomes. The population-specific differences in infant microbiome composition and function underscore the importance of studying microbiomes from people outside of wealthy, industrialized nations. Recognition of work on indigenous communities Research involving indigenous communities is needed for a variety of reasons including to ensure that scientific discoveries and understanding appropriately represent all populations and do not only benefit those living in industrialized nations. Special considerations must be made to ensure that this research is conducted ethically and in a non-exploitative manner. In this study we performed deep metagenomic sequencing on fecal samples that were collected from Hadza hunter-gatherers in 2013/2014 and were analyzed in previous publications using different methods ( 1, 2 ). A material transfer agreement with the National Institute for Medical Research in Tanzania ensures that stool samples collected are used solely for academic purposes, permission for the study was obtained from the National Institute of Medical Research (MR/53i 100/83, NIMR/HQ/R.8a/Vol.IX/1542) and the Tanzania Commission for Science and Technology, and verbal consent was obtained from the Hadza after the study’s intent and scope was described with the help of a translator. The publications that first described these samples included several scientists and Tanzanian field-guides as co-authors for the critical roles they played in sample collection, but as no new samples were collected in this study, only scientists who contributed to the analyses described here were included as co-authors in this publication. It is currently not possible for us to travel to Tanzania and present our results to the Hadza people, however we intend to do so once the conditions of the COVID-19 pandemic allow it.

Ultra-deep Sequencing of Hadza Hunter-Gatherers Recovers Vanishing Microbes (preprint)

The gut microbiome has been identified as a key to immune and metabolic health, especially in industrialized populations 1 . Non-industrialized individuals harbor more diverse microbiomes and distinct bacterial lineages 2 , but systemic under-sampling has hindered insight into the extent and functional consequences of these differences 3 . Here, we performed ultra-deep metagenomic sequencing and laboratory strain isolation on fecal samples from the Hadza, hunter-gatherers in Tanzania, and comparative populations in Nepal and California. We recover 94,971 total genomes of bacteria, archaea, bacteriophage, and eukaryotes, and find that 43% are novel upon aggregating with existing unified datasets 4,5 . Analysis of in situ growth rates, genetic pN/pS signatures, and high-resolution strain tracking reveal dynamics in the hunter-gatherer gut microbiome that are distinct from industrialized populations. Industrialized versus Hadza gut microbes are enriched in genes associated with oxidative stress, possibly a result of microbiome adaptation to inflammatory processes. We use phylogenomics to reveal that global spread of the spirochaete Treponema succinifaciens parallels historic human migration prior to its extinction in industrialized populations. When combined with a detailed definition of gut-resident strains that are vanishing in industrialized populations, our data demonstrate extensive perturbation in many facets of the gut microbiome brought on by the industrialized lifestyle. Recognition of work with indigenous communities Research involving indigenous communities is needed for a variety of reasons including to ensure that scientific discoveries and understanding appropriately represent all populations and do not only benefit those living in industrialized nations 3,6 . Special considerations must be made to ensure that this research is conducted ethically and in a non-exploitative manner. In this study we performed deep metagenomic sequencing on fecal samples that were collected from Hadza hunter-gatherers in 2013/2014 and were analyzed in previous publications using different methods 2,7 . A material transfer agreement with the National Institute for Medical Research in Tanzania specifies that stool samples collected are used solely for academic purposes, permission for the study was obtained from the National Institute of Medical Research (MR/53i 100/83, NIMR/HQ/R.8a/Vol.IX/1542) and the Tanzania Commission for Science and Technology, and verbal consent was obtained from the Hadza after the study’s intent and scope was described with the help of a translator. The publications that first described these samples included several scientists and Tanzanian and Nepali field-guides as co-authors for the critical roles they played in sample collection, but as no new samples were collected in this study, only scientists who contributed to the analyses described here were included as co-authors in this publication. It is currently not possible for us to travel to Tanzania and present our results to the Hadza people, however we intend to do so once the conditions of the COVID-19 pandemic allow it.

A 2-Gene Host Signature for Improved Accuracy of COVID-19 Diagnosis Agnostic to Viral Variants (preprint)

The continued emergence of SARS-CoV-2 variants is one of several factors that may cause false negative viral PCR test results. Such tests are also susceptible to false positive results due to trace contamination from high viral titer samples. Host immune response markers provide an orthogonal indication of infection that can mitigate these concerns when combined with direct viral detection. Here, we leverage nasopharyngeal swab RNA-seq data from patients with COVID-19, other viral acute respiratory illnesses and non-viral conditions (n=318) to develop support vector machine classifiers that rely on a parsimonious 2-gene host signature to predict COVID-19. Optimal classifiers achieve an area under the receiver operating characteristic curve (AUC) greater than 0.9 when evaluated on an independent RNA-seq cohort (n=553). We show that a classifier relying on a single interferon-stimulated gene, such as IFI6 or IFI44, measured in RT-qPCR assays (n=144) achieves AUC values as high as 0.88. Addition of a second gene, such as GBP5, significantly improves the specificity compared to other respiratory viruses. The performance of a clinically practical 2-gene RT-qPCR classifier is robust across common SARS-CoV-2 variants, including Omicron, and is unaffected by cross-contamination, demonstrating its utility for improving accuracy of COVID-19 diagnostics.

A 2-Gene Host Signature for Improved Accuracy of COVID-19 Diagnosis Agnostic to Viral Variants (preprint)

The continued emergence of SARS-CoV-2 variants is one of several factors that may cause false negative viral PCR test results. Such tests are also susceptible to false positive results due to trace contamination from high viral titer samples. Host immune response markers provide an orthogonal indication of infection that can mitigate these concerns when combined with direct viral detection. Here, we leverage nasopharyngeal swab RNA-seq data from patients with COVID-19, other viral acute respiratory illnesses and non-viral conditions (n=318) to develop support vector machine classifiers that rely on a parsimonious 2-gene host signature to predict COVID-19. Optimal classifiers achieve an area under the receiver operating characteristic curve (AUC) greater than 0.9 when evaluated on an independent RNA-seq cohort (n=553). We show that a classifier relying on a single interferon-stimulated gene, such as IFI6 or IFI44, measured in RT-qPCR assays (n=144) achieves AUC values as high as 0.88. Addition of a second gene, such as GBP5, significantly improves the specificity compared to other respiratory viruses. The performance of a clinically practical 2-gene RT-qPCR classifier is robust across common SARS-CoV-2 variants, including Omicron, and is unaffected by cross-contamination, demonstrating its utility for improving accuracy of COVID-19 diagnostics.

Using genomic epidemiology of SARS-CoV-2 to support contact tracing and public health surveillance in rural Humboldt County, California. (preprint)

During the COVID-19 pandemic within the United States, much of the responsibility for diagnostic testing and epidemiologic response has relied on the action of county-level departments of public health. Here we describe the integration of genomic surveillance into epidemiologic response within Humboldt County, a rural county in northwest California. Through a collaborative effort, 853 whole SARS-CoV-2 genomes were generated, representing [~]58% of the 1,449 SARS-CoV-2-positive cases detected in Humboldt County as of mid-March 2021. Phylogenetic analysis of these data was used to develop a comprehensive understanding of SARS-CoV-2 introductions to the county and to support contact tracing and epidemiologic investigations of all large outbreaks in the county. In the case of an outbreak on a commercial farm, viral genomic data were used to validate reported epidemiologic links and link additional cases within the community who did not report a farm exposure to the outbreak. During a separate outbreak within a skilled nursing facility, genomic surveillance data were used to rule out the putative index case, detect the emergence of an independent Spike:N501Y substitution, and verify that the outbreak had been brought under control. These use cases demonstrate how developing genomic surveillance capacity within local public health departments can support timely and responsive deployment of genomic epidemiology for surveillance and outbreak response based on local needs and priorities.

Upper airway gene expression reveals a more robust innate and adaptive immune response to SARS-CoV-2 in children compared with older adults (preprint)

Unlike other respiratory viruses, SARS-CoV-2 disproportionately causes severe disease in older adults and only rarely in children. To investigate whether differences in the upper airway immune response could contribute to this disparity, we compared nasopharyngeal gene expression in 83 children (<19-years-old; 38 with SARS-CoV-2, 11 with other respiratory viruses, 34 with no virus) and 154 adults (>40-years-old; 45 with SARS-CoV-2, 28 with other respiratory viruses, 81 with no virus). Expression of interferon-stimulated genes (ISGs) was robustly activated in both children and adults with SARS-CoV-2 compared to the respective non-viral groups, with only relatively subtle distinctions. Children, however, demonstrated markedly greater upregulation of pathways related to B cell and T cell activation and proinflammatory cytokine signaling, including TNF, IFN{gamma}, IL-2 and IL-4 production. Cell type deconvolution confirmed greater recruitment of B cells, and to a lesser degree macrophages, to the upper airway of children. Only children exhibited a decrease in proportions of ciliated cells, the primary target for SARS-CoV-2, upon infection with the virus. These findings demonstrate that children elicit a more robust innate and adaptive immune response to SARS-CoV-2 infection in the upper airway that likely contributes to their protection from severe disease in the lower airway.

Impaired immune signaling and changes in the lung microbiome precede secondary bacterial pneumonia in COVID-19 (preprint)

Secondary bacterial infections, including ventilator-associated pneumonia (VAP), lead to worse clinical outcomes and increased mortality following viral respiratory infections including in patients with coronavirus disease 2019 (COVID-19). Using a combination of tracheal aspirate bulk and single-cell RNA sequencing (scRNA-seq) we assessed lower respiratory tract immune responses and microbiome dynamics in 28 COVID-19 patients, 15 of whom developed VAP, and eight critically ill uninfected controls. Two days before VAP onset we observed a transcriptional signature of bacterial infection. Two weeks prior to VAP onset, following intubation, we observed a striking impairment in immune signaling in COVID-19 patients who developed VAP. Longitudinal metatranscriptomic analysis revealed disruption of lung microbiome community composition in patients with VAP, providing a connection between dysregulated immune signaling and outgrowth of opportunistic pathogens. These findings suggest that COVID-19 patients who develop VAP have impaired antibacterial immune defense detectable weeks before secondary infection onset.

Impaired antibacterial immune signaling and changes in the lung microbiome precede secondary bacterial pneumonia in COVID-19 (preprint)

Secondary bacterial infections, including ventilator associated pneumonia (VAP), lead to worse clinical outcomes and increased mortality following viral respiratory infections. Critically ill patients with coronavirus disease 2019 (COVID-19) face an elevated risk of VAP, although susceptibility varies widely. Because mechanisms underlying VAP predisposition remained unknown, we assessed lower respiratory tract host immune responses and microbiome dynamics in 36 patients, including 28 COVID-19 patients, 15 of whom developed VAP, and eight critically ill controls. We employed a combination of tracheal aspirate bulk and single cell RNA sequencing (scRNA-seq). Two days before VAP onset, a lower respiratory transcriptional signature of bacterial infection was observed, characterized by increased expression of neutrophil degranulation, toll-like receptor and cytokine signaling pathways. When assessed at an earlier time point following endotracheal intubation, more than two weeks prior to VAP onset, we observed a striking early impairment in antibacterial innate and adaptive immune signaling that markedly differed from COVID-19 patients who did not develop VAP. scRNA-seq further demonstrated suppressed immune signaling across monocytes/macrophages, neutrophils and T cells. While viral load did not differ at an early post-intubation timepoint, impaired SARS-CoV-2 clearance and persistent interferon signaling characterized the patients who later developed VAP. Longitudinal metatranscriptomic analysis revealed disruption of lung microbiome community composition in patients who developed VAP, providing a connection between dysregulated immune signaling and outgrowth of opportunistic pathogens. Together, these findings demonstrate that COVID-19 patients who develop VAP have impaired antibacterial immune defense weeks before secondary infection onset.

COVID-19 ARDS is characterized by a dysregulated host response that differs from cytokine storm and is modified by dexamethasone (preprint)

We performed comparative lower respiratory tract transcriptional profiling of 52 critically ill patients with the acute respiratory distress syndrome (ARDS) from COVID-19 or from other etiologies, as well as controls without ARDS. In contrast to a cytokine storm, we observed reduced proinflammatory gene expression in COVID-19 ARDS when compared to ARDS due to other causes. COVID-19 ARDS was characterized by a dysregulated host response with increased PTEN signaling and elevated expression of genes with non-canonical roles in inflammation and immunity that were predicted to be modulated by dexamethasone and granulocyte colony stimulating factor. Compared to ARDS due to other types of viral pneumonia, COVID-19 was characterized by impaired interferon-stimulated gene expression (ISG). We found that the relationship between SARS-CoV-2 viral load and expression of ISGs was decoupled in patients with COVID-19 ARDS when compared to patients with mild COVID-19. In summary, assessment of host gene expression in the lower airways of patients with COVID-19 ARDS did not demonstrate cytokine storm but instead revealed a unique and dysregulated host response predicted to be modified by dexamethasone.

Identification of a polymorphism in the N gene of SARS-CoV-2 that adversely impacts detection by a widely-used RT-PCR assay (preprint)

We identify a mutation in the N gene of SARS-CoV-2 that adversely affects annealing of a commonly used RT-PCR primer; epidemiologic evidence suggests the virus retains pathogenicity and competence for spread. This reinforces the importance of using multiple targets, preferably in at least 2 genes, for robust SARS-CoV-2 detection. Article Summary LineA SARS-CoV-2 variant that occurs worldwide and has spread in California significantly affects diagnostic sensitivity of an N gene assay, highlighting the need to employ multiple viral targets for detection.

The SARS-CoV-2 Spike mutation D614G increases entry fitness across a range of ACE2 levels, directly outcompetes the wild type, and is preferentially incorporated into trimers (preprint)

Early in the current pandemic, the D614G mutation arose in the Spike protein of SARS-CoV-2 and quickly became the dominant variant globally. Mounting evidence suggests D614G enhances viral entry. Here we use a direct competition assay with single-cycle viruses to show that D614G outcompetes the wildtype. We developed a cell line with inducible ACE2 expression to confirm that D614G more efficiently enters cells with ACE2 levels spanning the different primary cells targeted by SARS-CoV-2. Using a new assay for crosslinking and directly extracting Spike trimers from the pseudovirus surface, we found an increase in trimerization efficiency and viral incorporation of D614G protomers. Our findings suggest that D614G increases infection of cells expressing a wide range of ACE2, and informs the mechanism underlying enhanced entry. The tools developed here can be broadly applied to study other Spike variants and SARS-CoV-2 entry, to inform functional studies of viral evolution and vaccine development.

Epitopes targeted by T cells in convalescent COVID-19 patients (preprint)

There is growing evidence pointing to the protective role of T cells against COVID-19. Vaccines eliciting targeted T cell responses have the potential to provide robust, long-lasting immunity. However, their design requires knowledge of the SARS-CoV-2-specific epitopes that can elicit a T cell response and confer protection across a wide population. Here, we provide a unified description of emerging data of SARS-CoV-2 T cell epitopes compiled from results of 8 independent studies of convalescent COVID-19 patients. We describe features of these epitopes relevant for vaccine design, while indicating knowledge gaps that can, in part, be augmented using prior immunological data from SARS-CoV. The landscape of SARS-CoV-2 T cell epitopes that we describe can help guide SARS-CoV-2 vaccine development as well as future immunological studies. A web-based platform has also been developed to complement these efforts.

AI aided design of epitope-based vaccine for the induction of cellular immune responses against SARS-CoV-2 (preprint)

The heavy burden imposed by the COVID-19 pandemic on our society triggered the race towards the development of therapies or preventive strategies. Among these, antibodies and vaccines are particularly attractive because of their high specificity, low probability of drug-drug interaction, and potentially long-standing protective effects. While the threat at hand justifies the pace of research, the implementation of therapeutic strategies cannot be exempted from safety considerations. There are several potential adverse events reported after the vaccination or antibody therapy, but two are of utmost importance: antibody-dependent enhancement (ADE) and cytokine storm syndrome (CSS). On the other hand, the depletion or exhaustion of T-cells has been reported to be associated with worse prognosis in COVID-19 patients. This observation suggests a potential role of vaccines eliciting cellular immunity, which might simultaneously limit the risk of ADE and CSS. Such risk was proposed to be associated with FcR-induced activation of proinflammatory macrophages (M1) by Fu et al. 2020 and Iwasaki et al. 2020. All aspects of the newly developed vaccine (including the route of administration, delivery system, and adjuvant selection) may affect its effectiveness and safety. In this work we use a novel in silico approach (based on AI and bioinformatics methods) developed to support the design of epitope-based vaccines. We evaluated the capabilities of our method for predicting the immunogenicity of epitopes. Next, the results of our approach were compared with other vaccine-design strategies reported in the literature. The risk of immuno-toxicity was also assessed. The analysis of epitope conservation among other Coronaviridae was carried out in order to facilitate the selection of peptides shared across different SARS-CoV-2 strains and which might be conserved in emerging zootic coronavirus strains. Finally, the potential applicability of the selected epitopes for the development of a vaccine eliciting cellular immunity for COVID-19 was discussed, highlighting the benefits and challenges of such an approach.

Upper airway gene expression differentiates COVID-19 from other acute respiratory illnesses and reveals suppression of innate immune responses by SARS-CoV-2 (preprint)

We studied the host transcriptional response to SARS-CoV-2 by performing metagenomic sequencing of upper airway samples in 238 patients with COVID-19, other viral or non-viral acute respiratory illnesses (ARIs). Compared to other viral ARIs, COVID-19 was characterized by a diminished innate immune response, with reduced expression of genes involved in toll-like receptor and interleukin signaling, chemokine binding, neutrophil degranulation and interactions with lymphoid cells. Patients with COVID-19 also exhibited significantly reduced proportions of neutrophils and macrophages, and increased proportions of goblet, dendritic and B-cells, compared to other viral ARIs. Using machine learning, we built 26-, 10- and 3-gene classifiers that differentiated COVID-19 from other acute respiratory illnesses with AUCs of 0.980, 0.950 and 0.871, respectively. Classifier performance was stable at low viral loads, suggesting utility in settings where direct detection of viral nucleic acid may be unsuccessful. Taken together, our results illuminate unique aspects of the host transcriptional response to SARS-CoV-2 in comparison to other respiratory viruses and demonstrate the feasibility of COVID-19 diagnostics based on patient gene expression.

Clinical features, diagnostics, and outcomes of patients presenting with acute respiratory illness: a comparison of patients with and without COVID-19 (preprint)

BackgroundEmerging data on the clinical presentation, diagnostics, and outcomes of patients with COVID-19 have largely been presented as case series. Few studies have compared these clinical features and outcomes of COVID-19 to other acute respiratory illnesses. MethodsWe examined all patients presenting to an emergency department in San Francisco, California between February 3 and March 31, 2020 with an acute respiratory illness who were tested for SARS-CoV-2. We determined COVID-19 status by PCR and metagenomic next generation sequencing (mNGS). We compared demographics, comorbidities, symptoms, vital signs, and laboratory results including viral diagnostics using PCR and mNGS. Among those hospitalized, we determined differences in treatment (antibiotics, antivirals, respiratory support) and outcomes (ICU admission, ICU interventions, acute respiratory distress syndrome, cardiac injury). FindingsIn a cohort of 316 patients, 33 (10%) tested positive for SARS-CoV-2; 31 patients, all without COVID-19, tested positive for another respiratory virus (16%). Among patients with additional viral testing, no co-infections with SARS-CoV-2 were identified by PCR or mNGS. Patients with COVID-19 reported longer symptoms duration (median 7 vs. 3 days) and were more likely to report fever (82% vs. 44%) fatigue (85% vs. 50%) and myalgias (61% vs 27%); p<0.001 for all comparisons. Lymphopenia (55% vs 34%, p=0.018) and bilateral opacities on initial chest radiograph (55% vs. 24%, p=0.001) were more common in patients with COVID-19. Patients with COVID-19 were more often hospitalized (79% vs. 56%, p=0.014). Of 186 hospitalized patients, patients with COVID-19 had longer hospitalizations (median 10.7d vs. 4.7d, p<0.001) and were more likely to develop ARDS (23% vs. 3%, p<0.001). Most comorbidities, home medications, signs and symptoms, vital signs, laboratory results, treatment, and outcomes did not differ by COVID-19 status. InterpretationWhile we found differences in clinical features of COVID-19 compared to other acute respiratory illnesses, there was significant overlap in presentation and comorbidities. Patients with COVID-19 were more likely to be admitted to the hospital, have longer hospitalizations and develop ARDS, and were unlikely to have co-existent viral infections. These findings enhance understanding of the clinical characteristics of COVID-19 in comparison to other acute respiratory illnesses.